Experimental estimation of entanglement at the quantum limit 
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Entanglement is the central resource of quantum information processing and the precise charac- 
terization of entangled states is a crucial issue for the development of quantum technologies. This 
leads to the necessity of a precise, experimental feasible measure of entanglement. Nevertheless, such 
measurements are limited both from experimental uncertainties and intrinsic quantum bounds. Here 
we present an experiment where the amount of entanglement of a family of two-qubit mixed photon 
states is estimated with the ultimate precision allowed by quantum mechanics. 
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Entanglement is the central resource of quantum in- 
formation processing and the precise characterization of 
entangled states is a crucial issue for the development of 
quantum technologies. In turn, quantification and detec- 
tion of entanglement have been extensively investigated, 
see 0, Q, 0] for a review, and different approaches have 
been developed to extract the amount of entanglement of 
a state from a given set of measurement results [i],0,0,0]- 
Of course, in order to evaluate the entanglement of a 
quantum state one may resort to full quantum state to- 
mography Q which, however, becomes impractical in 
higher dimensions and may be affected by large uncer- 
tainty [9(| . Other methods, requiring a reduced number 
of observables, are based on visibility measurements II di . 
Bell' tests 11, 12, 13 1, entanglement witnesses [I3.ll5lllq 



or are related to Schmidt number 
found an experimental application 



in the presence of decoherence effects 
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lJUMJp, [H, also 

Any quantitative measure of entanglement corresponds 
to a nonlinear function of the density operator and thus 
cannot be associated to a quantum observable. As a con- 
sequence, ultimate bounds to the precision of entangle- 
ment measurements cannot be inferred from uncertainty 
relations. Any procedure aimed to evaluate the amount 
of entanglement of a quantum state is ultimately a pa- 
rameter estimation problem, where the value of entan- 
glement is indirectly inferred from the measurement of 
one or more proper observables [25!]. An optimization 
problem thus naturally arises, which may be properly 
addressed in the framework of quantum estimation the- 
ory 0, [2?J , which provides analytical tools to find the 
optimal measurement and to derive ultimate bounds to 
the precision of entanglement estimation. 

Suppose one has a family of quantum states g e la- 
beled by the value of entanglement, say negativity, and 



wants to estimate e from the outcomes x — ( x ii --^n) 
of N repeated measurements of the (generalized) observ- 
able described by a Positive Operator Valued Measure- 
ment (POVM) II X , J2x^x = Am/ inference strat- 
egy amounts to find an estimator, i.e. a map e(x) 
from the experimental sample to the parameter space. 
According to the Cramer-Rao theorem the precision of 
any estimation procedure, i.e. the variance of any un- 
biased estimator based on the measurement of is 
bounded by the inequality Var(e) > [JVFJ -1 , where 
F e — J2 x P( x \ e ) ft* lnp(a;|e)] 2 is the Fisher information 
and p(x\e) = Tr[g e H x ] is the conditional probability of 
getting the outcome x when the actual value of entangle- 
ment is e. Upon maximizing the Fisher information over 
all the possible quantum measurements we arrive at the 
quantum Fisher information (QFI) F e < H t — Tr[X 2 p e ] 
expressed in terms of the symmetric logarithmic deriva- 
tive L e , i.e. the selfadjoint operator satisfying the equa- 
tion d e p e = ^(L e p > : + p t L t ). The ultimate bounds to 
precision are thus determined by the quantum Cramer- 
Rao bound (QCRB) Var(e) > [NF^" 1 > [NHe]- 1 . The 
meaning of QCRB is that quantum mechanics does not 
allow entanglement estimation with arbitrary precision. 
In turn, QCRB represents the ultimate bound to the pre- 
cision of any procedure aimed to estimate the amount of 
entanglement of a state belonging to the family g e . In 
order to optimally estimate entanglement we thus need 
i) a measurement with Fisher information F e = H e equal 
to the QFI and ii) an estimator saturating the Cramer- 
Rao bound [2^|. In [25}, bounds to precision have been 
evaluated for several classes of pure and mixed quan- 
tum states. Here we demonstrate experimentally for the 
first time that optimized correlations measurements al- 
lows for the estimation of entanglement with the ultimate 
precision imposed by quantum mechanics. In particular, 
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we present the results of an experiment to estimate the 
amount of entanglement (negativity) of two-qubit pho- 
ton states. This represents a substantial advance, paving 
the way for further progresses. In fact, with a judicious 
choice of correlation measurements one can devise a pro- 
cedure to optimally estimate entanglement for a generic 
class of two-photon entangled states. 

The family of entangled states we are dealing with is 
made of polarization entangled photon pairs obtained by 
coherently superimposed type-I parametric downconver- 
sion (PDC) generated in two BBO crystals [lj|. The 
experimental set-up is schematically depicted in Fig. 
[TJ A CW Argon pump laser beam with wavelength 
A = 351nm is filtered with a dispersion prism and then 
passes through a Glan-Thompson prism with horizontal 
axis of transmission that selects the polarization \H). An 
halfwave plate WPO rotates the polarization by the angle 
</>, which in turn determines the amount of entanglement 
in the output state. PDC light is generated by two thin 
type-I BBO crystals (l=lmm), positioned with the planes 
that contain optical axes orthogonal to each other. PDC 
occurs only in crystal 1(2) if the polarization of the pump 
beam is horizontal(vertical). 
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FIG. 1: (Color online) Experimental setup to generate po- 
larization entangled photon pairs with variable entanglement 
and estimate its value with the ultimate precision allowed by 
quantum mechanics. 

The crystals are cut for collinear frequency degener- 
ate phase-matching at working wavelength and the phase 
shifts due to ordinary and extraordinary path in the 
crystals are compensated by rotating the quartz plates 
QP. To maintain stable phase-matching conditions, BBO 
crystals and QP are placed in a closed box which is kept 
heated at fixed temperature. Overall, the output states 
are described by the family of density matrices g e = 
pl^/vKVV'l + (1 ~P)D<p which represent the mixing of the 
ideal pure entangled states |Vv>) = cos <j)\HH)+sva.<j)\VV) 
(negativity e = sin 20) with a small fraction of a separa- 
ble mixture D^> = cos 2 (j>\HH){HH\ + sin 2 cj>\VV){VV\. 
The angle may be tuned at will upon rotating the 
waveplate WPO, whereas the mixing parameter p comes 
from the decoherence mechanisms occurring in the exper- 
imental setup. These are mostly due to fluctuations of 
the relative phase between the two polarization compo- 
nents, which themselves derive from residual temperature 
fluctuations: in order to obtain different mixed states 



these fluctuation are eventually magnified by eliminat- 
ing the isolating box. The entanglement (negativity) of 
the state g e is given by e — p sin 20 and its purity by 
fi = 1 — (1 — p 2 ) sin 20. Upon inverting these relations 
and express the family of states in terms of e and /i or 
p we may evaluate the QFI for entanglement estimation, 
which turns out to be a function of e only, H e = (1— e 2 ) -1 . 
In the following, we describe our detection strategy and 
show it allows entanglement estimation with precision 
saturating the QCRB independently on the purity. 

After the crystal, the pump is stopped by a filter 
(UVF), and the biphoton field is split on a non-polarizing 
50-50 beam splitter (BS). The measurement stage con- 
sists in projecting the beams on vertical polarizers af- 
ter passing through halfwave plates (WP1, WP2). After 
spectral selection by an interference filters (IF) centered 
at the degeneracy 702nm (FWHM=3nm) , biphotons are 
collimated by a short focal length lens and detected on 
SPAD (Dl, D2). Electrical signals from detectors are 
registered with a coincidence scheme (CC) with a time 
window of coincidences set to Ins. Overall, the measure- 
ment scheme is described by projection measurements 
onto two-qubit states 

n t (a,/3) = \a + s^){a + s~\ ® \p + s'^)(p + s'^\ 

where t = {s + 2s'}, s, s' = 0, 1. The polarization angles 
a, (3 are set by Glen- Thompson polarizers, whereas the 
rotations of ±7r/2 are obtained by putting the half-wave 
plates WP1 and WP2 at ±22,5° (mounted on precision 
rotation stages with high resolution and fully motor con- 
trolled). 

Let us first illustrate precision analysis assuming the 
generation of the pure states |Vv>)- In this case neg- 
ativity is given by e = sin 20 and the estimation of 
entanglement reduces to a measurement of coincidence 
rates in a two-particle interferometer setting 0. In- 
deed, upon inspecting the expression of the probabili- 
ties p t (e;a,f3) = (^|n t (a,/3)|Vty), t = 0,1,2,3 one finds 
out that an unbiased estimator for the negativity can 
be written as e = V(a, (3) esc 2a esc 2(3 — cot(2a) cot(2/3), 
where V(a,/3) = po — p\ — P2 + P3 is the expected 
value of two-qubit quantum correlations (QC). The cor- 
responding Fisher information is thus given by F e = 
J2t Pt{ e 'i a , P){dt logpt(e; a, /?)] 2 and it equals the QFI, 
H e , for a — ±7r/4, (3 = ±7r/4. In other words e = 
V r (±7r/4, ±7r/4), which can be measured with the exper- 
imental setup described above, are good candidates for 
being optimal estimators of entanglement. In practice, 
at fixed 0, a and (3 and for each measurement run j = 
1, .., M one records the vector kj = {ko t j, fcij, k2j, fcaj}, 
where k t j = ktj(a,/3), is the number of coincidence 
counts, as measured for the given set of parameters by 
the coincidence circuit during a single time window (10 
seconds). For large values of the total number of coin- 
cidences Kj = ^2 t k t .j (determined irrespectively of the 



3 



polarizers' orientation in the j — th run), the expected 
value of the coincidence rate k t ,j(a, (3) / Kj converges to 
the probability Pt(e; a, (3) and the estimator can thus be 
written in terms of the coincidences' vector e = e(kj). In 
our implementation we have performed M — 30 measure- 
ments of the coincidence vector for fixed values a = — 7r/4 
and (3 = tt/4. 

For finite KjS the uncertainty in the estimation of en- 
tanglement are mostly due to fluctuations 5k t in the co- 
incidence counts kt,j around their average values (kt) = 
y"V k t j/M. Using standard error propagation with the 
derivatives d t = d/dk t evaluated for k t = (kt), and 
assuming independence among fluctuations at different 
angles, we have Var(e) = J2t \dt£\ 2 8k 2 = 4[((fc ) + 
(k 3 )) 2 (5k 2 + 5k 2 ) + ((k 1 ) + (k 2 )) 2 (Sk 2 + 6k 2 )]/(Kf. If we 
now assume that the counting processes have a Poisso- 
nian statistics, i.e. 8k 2 = (k t ) 2 , then it is straightforward 
to prove that 

Var(e) = 4(fc + fa)(fei + fa)/ (K) 3 = (1 - e 2 )/ (K) 

i.e. QC measurements allow for optimal estimation of 
entanglement with precision at the quantum limit. Since 
the QFI is given by H e = (1 — e 2 ) -1 for a wide range of 
two-qubit families of states [25|, the above calculations 
suggest that this is a general result. In other words, 
given a source emitting polarization two-qubit states with 
coincidence counting statistics satisfying the Poissonian 
hypothesis, then the experimental setup of Fig [T] allows 
for optimal estimation of entanglement at the quantum 
limit by means of a QC estimator. 

This is indeed the case for the experimental setup of 
Fig. [1] and the family of output states g e where, upon 
evaluating the probabilities pt(e;a,j3) — Tr[g e 11* (a,/?)], 
one sees that e = V(— tt/4, 7r/4) is still an optimal (un- 
biased) estimator of entanglement. In practice, we have 
collected M = 30 repeated acquisitions of coincidence 
vector kj — {fay, fey, faj, faj}, then we have random- 
ized the composition of kj over the sequence of mea- 
surements to avoid spurious correlations, and finally we 
have estimated entanglement as the sample mean (e) = 
e(kj)/M. The corresponding uncertainty has been 
evaluated by the sample variance Var(e) = YljiK^j) ~ 
(e)] 2 /(M - 1). In order to compare the estimated value 
of entanglement with the actual one we need to esti- 
mate also the additional parameter p, quantifying the 
amount of mixing introduced by decoherence processes. 
An unbiased estimator p for this parameter may be ob- 
tained by measuring QC with a different set of angles, 
e.g. upon collecting the coincidences rj = kj(a = (3 = 0) 
to form p(rj,kj) = ~e(kj)Rj / 'yi^J(i — r 3 j), where 
Rj = ^2 t r ti j is the total number of coincidences with 
the four orientations a,j3 = 0, 7r/2. The actual ("true") 
value of negativity is then inferred as = (p) sin 2</>, i.e. 
using the knowledge of the rotation angle of the wave- 
plate WPO and the estimation of the mixing parameter. 




FIG. 2: Estimation of entanglement at the quantum limit. 
The plot shows the estimated value of entanglement (e) as a 
function of the actual one et. The uncertainty bars on (e) 
denotes the quantity ^Vax(e) X (K), i.e. the square root 
of the sample variance multiplied by the average number of 
total coincidences (K). The gray area corresponds to val- 

— 1/2 

ues within the inverse of the Fisher information et ± H tt 
Uncertainty bars on the abscissae correspond to fluctuations 
Se t = \/Var(p) sin 2(f) in the determination of e t , due to fluc- 
tuations in the estimation of the mixing parameter. 

In Fig. [5] we show the estimated value of en- 
tanglement as a function of the actual one for seven 
values of <f> = 10°, 15°, 20°, 28°, 40°, 45°, 45° of the 
WPO rotation angle (corresponding to estimated mix- 
ing (p) = 0.85,0.88,0.88,0.85,0.92,0.93,0.97 respec- 
tively). The uncertainty bars on (e) denotes the quan- 
tity y/Vai(e) x (K), i.e. the square root of the sample 
variance multiplied by the average number of total co- 
incidences (K), This is in order to allow a direct com- 
parison with the Cramer-Rao bound in term of the in- 
verse of the Fisher information (the gray area). Uncer- 
tainty bars on the abscissae correspond to fluctuations 
Set = -\/Var(p) sin 2<f> in the determination of et , due to 
fluctuations in the estimation of the mixing parameter. 
As it is apparent from the plot entanglement is estimated 
with precision at the quantum limit for any value of the 
the rotation angle <f>. Notice that this conclusion is robust 
against the fact that the statistics is not exactly Poisso- 
nian: in the left panel of Fig. [3] we show the Fano factor 
for the four recorded coincidence counts kj and the six 
values of <j) employed in the experiment. 

Our statistical model g e may be checked for consis- 
tency on the basis of the recorded data themselves, and 
other possible models to describe decoherence of our fam- 
ily of states are ruled out as they cannot fit the experi- 
mental sample. As for example, if one tries to describe 
the output from our source by the family of (depolar- 
ized) Werner states g' e = p\ip</>) (4><j>\ + \(l — p)& ® 1, than 
one sees from the expression of the coincidence proba- 
bility p' t (e;a,P) = Tr[g' e H t (a,f3)] that unbiased estima- 
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tors for the mixing parameters and the negativity may 
be expressed as p' = 7(0,0), e' = -\ + ±P(0,0) + 
V(— it/4, 7r/4). These may be written in terms of the co- 
incidence vectors k and r as p'{vj) — (r rQ j — r^j — r 2 j + 
r 3 ,j)/Rj and e'ir^kj) = -± + ±p'(r 3 ) + (fc ,j - h d - 
&2J +ks t j)/Kj. Upon evaluating the corresponding sam- 
ple means and variances one realizes that the model is 
incompatible with the observed data. This is illustrated 
in the right panel of Fig|3] where we report the estimated 
value of entanglement as a function of the actual one 
assuming, for the description of the output signals, the 
families g e (top plot) and g' e (bottom plot). Here the 
uncertainty bars denote the 3c confidence interval and 
thus it is apparent that g' t cannot fit the data. We also 
performed measurements of the coincidence counts for a 
complete set of polarization angles and check the model 
g t by full two-qubit state tomography 
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FIG. 3: (Color online) Left: Fano factors of the coincidence 
counts kj, j — 0, 1, 2, 3. Each groups contains the Fano factor 
for the seven values of <j> reported in the text. Right: esti- 
mated value of entanglement as a function of the actual one 
assuming, for the description of the output signals, the fami- 
lies g t (top plot) and g' e (bottom plot). The uncertainty bars 
stays for to the 3<r confidence interval. 

In conclusion, we have suggested and demonstrated a 
measurement scheme based on quantum correlation mea- 
surements to optimally estimate entanglement for a fam- 
ily of two-photon entangled states. Our procedure is self- 
consistent and allows to estimate the amount of entan- 
glement with the ultimate precision imposed by quantum 
mechanics. With an appropriate choice of correlation 
measurements our results may be extended to a generic 
class of two-photon entangled states. The statistical re- 
liability of our method suggests a wider use in precise 
monitoring of external parameters assisted by entangle- 
ment. 
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